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Abstract 

In this article, the problem of the charged harmonic plus an inverse harmonic 
oscillator with time-dependent mass and frequency in a time-dependent elec- 
tromagnetic field is investigated. It is reduced to the problem of the inverse 
harmonic oscillator with time-independent parameters and the exact wave 
function is obtained. 
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I. INTRODUCTION 



A time-dependent quantum system with a Hamiltonian which contains a singular term 

P 2 m (t) 9/ x „2 C 

is one of the rare examples admitting exact solution and has been subject of great interest in 
a variety of contexts [1-10], because of its various applications in different areas of physics. 
For example, it is used for constructing exactly solvable models of interacting many-body 
systems [1,2], for modelling polyatomic-molecules [3] and for coherent-state formalism [11]. 
The Hamiltonian in Eq. (1) has been studied and the exact solution was obtained by 
using the Lewis and Riesenfeld invariant method and by making canonical and unitary 
transformations [5-9]. Although the exact solution of the time-dependent charged harmonic 
oscillator in a time- varying electromagnetic field is known [12], the charged harmonic plus 
inverse harmonic oscillator problem in a time-varying electromagnetic field has not yet been 
solved exactly. 

In this letter we take the following Hamiltonian into account 



H 



p-gA(t)l m (t) C 7/ , .0 



2m(t) 2 w ^ m(t)q 2 ^ v ' dt 

where we have set h — 1 = c. By making some transformations, Eq. (2) is reduced to 

£ + £ = </<oi, (3) 

where f(t) is a time- dependent function. Finally, Eq. (3) is integrated. 



II. METHOD 

The Hamiltonian for the charged harmonic plus inverse harmonic oscillator constrained 
to move in the x — y plane is given 



H 



i 2 

p-qA(t) m (t) 2 . w , 2 ^ C lus .d /A , 



2m(t) 2 m(t)(x 2 + y 2 ) dt 



d , d . 

where C is any constant and p = — — ei — i— e 2 . 

ot ay 

In this letter, we restrict ourselves to the case where the electromagnetic potentials are 
given by 

A = ^-(ye 1 -xe 2 ), 0(f) = 0, (5) 

where (ei, e 2 , e 3 ) are the unit vectors of the cartesian basis. With this choice, the magnetic 
and electric fields turn out to be 

B = B(t)e 3 , (6) 
1 dB(t) ? 

where (p, 0, e 3 ) form the cylindrical basis. Hereafter, we will denote cu(t), m(t) and 5(f) by 
ou, m and 5, respectively, for brevity. By substituting Eq. (5) into Eq. (4), the Hamiltonian 
becomes 

Px 2 +P y 2 qB m( 2 q 2 B 2 \ 2 2 C .9 

Now we want to make a transformation to get rid of the cross-term in Eq. (8). To do this, 
an orthogonal transformation is introduced [12]. 

x\ / cos (3(t) -sin /3(f) \ fx' 
y) \sm(3(t) cos /3(f) ) \y' 
The time- dependent transformation (9) implies the following transformations for the 
quadratic terms 

Px +Py ~>Px' +Py> , x +y ^x +y , (10) 

and for the cross-term and time derivative 

ypx - xp y -»• ypx' - x'py, (11) 

§- t ^§- t +^(t)[y'px'-xpy'}. (12) 



(9) 



To eliminate the cross-term from the Eq. (8), (3 is chosen as 

qB 



4m' 



(13) 



where dot denotes time derivation as usual. In this case, Eq. (8) becomes 

q 2 B 2 \ ^ C .0 



2m 2 \ Am 2 ) { U f m(x' 2 + y' 2 ) ' dt' 



(14) 



Eq. (14) does not depend on the cross-term and the square of the angular frequency is 

shifted by an amount of q 2 B 2 /(Am 2 ). 

To remove the harmonic part, we make the transformation 



pj -> A = Px + ia(t)x', py -> P 2 = Py + ia(t)y', 



(15) 



(16) 
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where a(t) is a time-dependent function to be determined later. Since p = — i— ei — i— e 2 , 

ox ay 

we observe that the transformation (15) is equivalent to the multiplication of the wave 
function by exp -a(t)(x' 2 + y' 2 )/2 . 

Under the transformation (15), the time derivative operator becomes 



d d ( Qi 2 + Q2 2 



Substituting Eqs. (15), (16) and (17) in (14), we get 

C 



(17) 



where 



Qi +Q2 2 



— — + i^(QiA + O a P 2 )=<^ + ^, (18) 



m(Q 1 +Q 2 ) m 



dt m 



„. . m / 2 q 2 B 2 \ a 2 (t) .a(t) 

n 2 (t) = — \uj 2 + - + i-^-. 

w 2 \ 4m 2 J 2m 2 

The time-dependent functions a(t) should be chosen so that harmonic term in the Hamil- 
tonian vanishes. 



m q 2 ^ a 2 ^ a^ = 

2 [UJ + 4m2 j 2m + * 2 U ' 



(19) 



Eq. (18) is simplified after the elimination of the harmonic part. 

~ 2 - 2 

Pi +P 2 c 



2m m [Qi+Q* 



+ ia(t) [4A + 4A 



d a(t) 
at m 



(20) 



Now a velocity-dependent interaction term is contained in Eq. (20). To get rid of this term, 
we need another transformation. 



Pi — > Pi = v(t)Pi, P2 -> P2 = ti(t)P 2 



(21) 



Qi ~* Qi — ' ~ * 



tity 

then the time derivative operator becomes 



(22) 



(23) 



If /x(t) is chosen so that 



a(t) 



AW 



(24) 



then the velocity dependent interaction term vanishes. The new equation is just the equation 
of a particle under an inverse harmonic potential and with a time-dependent mass. 



P x +P 2 + — — = i2mfi 2 (t)— + 2/i 2 (t)a(t) = k 2 . 

Ql + <?2 



(25) 



Since left-hand side and right-hand side of (25) depend on different parameters, the two 
sides are equal to a constant, k 2 . We can rewrite Eq. (25) as a wave equation rather than 
an operator equation 



d 2 ^(x,y,t) d 2 ^(x,y,t) 2C 



+ 



ty(x : y, t) = i2mfx 



,dV(x,y,t) 



dQ'i 2 dQ' 2 2 ' Q? + gf dt 

Using the separation of variables technique, we can solve Eq. (26). 
Let us write * = R(Q[, Q' 2 )T(t) 

rt k 2 + 2^{t')a{t') J 



+ 2fi 2 a^(x,y,t). (26) 



T(t) = exp 
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Jo 



2m(t')fi 2 (t>) 



-dt' 



(27) 



Coordinate-dependent part of Eq. (26) can be decomposed in polar coordinates 

Q 2 ^ + + (* 2 Q 2 - 2C - n 2 )Z = (28) 

§ - - 2 e, (29) 

where Q 2 = Qf + Q' 2 2 , Q\ = Q cos 9, Q' 2 = Qsin9 and R(Q[, Q' 2 ) = Z(Q)&(6). 
The relation between 9 and x, y is given 

tanfl = f = cosftr + sin/fr 
— sin + cos 

where /3 was defined in Eq. (13). 

If we introduce v 2 = 2C + n 2 , we see that Eq. (28) is just the Bessel differential equation 
of order v [13]. The solutions of (28,29) 

Z(Q) = AJ u (kQ) + BN u (kQ) (31) 
e = Cexp(±in0), (32) 

where J u (kQ'), N u {kQ') are the Bessel function of the first and the second kind respectively. 
The constants A, B and C are determined from the boundary conditions. Note that v and 
n are not the same. If the inverse harmonic potential term goes to zero, they become equal 
to each other. 

Finally, the complete exact solution of Eq. (4) is given 



m = 



exp[a(t)(x 2 + y 2 ) ±in9-if] (33) 



r l k 2 + 2u 2 (t')a(t') 

where f(t) = / " v ' v ' dt' and a, /i, 9 were defined Eqs. (19), (24) and (30) 

respectively. 

In summary, the exact solution of the problem of the charged harmonic plus an inverse 
harmonic oscillator with time-dependent mass and frequency in a time- dependent electro- 
magnetic field is obtained by making some appropriate transformation on momentum and 
coordinate operators. 
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